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Abstract 

In this paper we construct polynomial lattice rules which have, in some sense, 
small gain coefficients using a component-by-component approach. The gain co- 
efficients, as introduced by Owen, indicate to what degree the method improves 
upon Monte Carlo. We show that the variance of an estimator based on a scram- 
bled polynomial lattice rule constructed component-by-component decays at a rate 
of Ar-(2«+i)+5^ for all -5 > 0, assuming that the function under consideration has 
bounded variation of order a and where denotes the number of quadrature points. 
An analogous result is obtained for Korobov polynomial lattice rules. It is also 
established that these rules are almost optimal for the function space considered 
in this paper. Furthermore, we discuss the implementation of the component-by- 
component approach and show how to reduce the computational cost associated 
with it. Finally, we present numerical results comparing scrambled polynomial lat- 
tice rules and scrambled digital nets. 
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1 Introduction 

Quasi-Monte Carlo rules '^n=i fi^n), Xi, . . . ,Xn G [0, l]**, are equal weight integration 
formulas used to approximate integrals over the unit cube [0, 1]**, where s is typically 
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large. One can roughly divide quasi- Monte Carlo rules into lattice rules, see e.g. [T81I27] . 
and digital nets, see e.g. [3 HHj. In this paper we focus on digital nets, the construction 
of which is based on linear algebra over finite fields, see [71 [18]. In particular, we are 
interested in a special case of digital nets, so-called polynomial lattice rules, which are 
constructed using polynomials over finite fields; polynomial lattice rules were introduced 
in dg, see also [SUZIOK]. 

Studying the approximation of integrals using quasi-Monte Carlo methods, one wants 
to have information on the resulting integration errors. However, depending on the inte- 
grand under consideration, estimates of integration errors might be very conservative or 
unknown; a possible remedy to this problem is randomization, which allows us to obtain 
statistical information on integration errors, [23]. Popular choices of randomization meth- 
ods are digital shifts, see e.g. [HI E], and scrambling as introduced by Owen [23], see also 
[3 ini [HI [211 [2S1 [301 ISIl [321 [33] • in this paper, we focus on scrambling. In particular, we 
are interested in the variance of the estimator 

h=0 

where the points {yh}h=o obtained by applying the scrambling algorithm to a poly- 
nomial lattice rule. Notice that /(/) is an unbiased estimator of f{x)dx, that is, 

EW)) = /[o,i]» /Nd^' see [23]- 

The variance of the estimator given in Equation admits the representation, see 

Var(J(/)) = l ^^^i(f^^ (2) 

«eNg\{0} 

where N is the number of quadrature points. Equation (|2]) holds for any estimator 
obtained by applying the scrambling algorithm to a point set {x^}^^^^ such that x^ G 
[0, 1)'^. Here, the values Ti are the so-called gain coefficients which depend only on the 
quadrature points and the values cri{f) depend only on the integrand /. They are derived 
from the crossed and nested Anova decomposition of /, see [21], and can be expressed in 
terms of Haar coefficients of the function /, see [21], or also as a sum of certain Walsh 
coefficients of /, see [3 Section 13.2]. In this sense. Equation (|2]) shows that Var(J(/)) 
can be expressed as a weighted sum of gain coefficients, where we interpret the <yf{f) as 
weights. 

In our investigations we consider a space of functions for which o"j(/) has a certain 



3 



rate of decay. More precisely, for < a < 1 we introduce a norm of the form 

||/|U = sup6"l'lv,(/), (3) 

where \l\i = h + ■ — \- h for Z = (/i, . . . , Z^). We show that a is related to the smoothness 
of / in the following sense: 

If / G L2{[0, lY) has bounded variation of order a, then ||/||q, < oo. 

See Corollary 12.11 for details. 

From ([3]) we obtain ai{f) < and by substituting this formula into we 
obtain 

Var(/(/))<^ Yl r.&-^""Hl/lll (4) 

«GNg\{0} 

To construct polynomial lattice rules of high quality we use 

«eNg\{o} 

as quality criterion. Notice that the sum ([5]) only depends on the quadrature points and 
not on the function /. We show that (jS]) has a simple closed form for any < a < 1 
which can easily be computed if the quadrature points are a digital net. The case a = 
needs to be excluded since in this case ([5]) is infinite. 

Our aim is to find polynomial lattice rules for which the weighted sum of gain coef- 
ficients ([5]) is minimized. It is known from [2U [25l [31] that a small quality parameter 
t of a digital (t, m,s)-net yields small gain coefficients. In fact, one has = for all 
Z G Nq \ {0} such that |Z|i < m — t. Here, on the other hand, we aim to minimize (I5j) 
since it can be used to bound the variance of the estimator Var(J(/)). In other words, 
we minimize the upper bound on the variance (jlj) for all functions / with ||/||q, < oo over 
the class of polynomial lattice rules. 

We introduce additional parameters 7 = {'Jj)i<j<s in the norm in the sense of 
[25] . In this case the criterion depends on the additional parameters 7 = {'Jj)i<j<s, 
in which case a small quality parameter t does not necessarily yield the smallest possible 
gain coefficients anymore. In such a situation component-by-component constructions [28] 
have proven useful since one can then optimize the quadrature points also with respect to 
the 7j. This is also the approach taken here. More precisely, we show that by constructing 
polynomial lattice rules component-by-component one obtains a convergence of 



Var(/(/)) = 0(iV-(2"+^)+'') for any 6 > 0. 
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Apart from S, which can be arbitrarily close to 0, this rate is best possible as shown in 
Section [6] for a large class of randomized algorithms. Further, if YlJLilj < oo, then the 
bound (jl]) does not depend on the dimension s. This result is stronger than what can be 
obtained for {t, m, s)-nets, since the increase of the t value of the form t ^ s prevents one 
from obtaining a bound independent of the dimension only assuming that Yl'jLilj < 
This condition is also necessary, see [52] for a result on a related space. Hence the rules 
we construct here are simultaneously optimal in terms of the convergence rate as well as 
in terms of their dependence on the dimension. 

Notice that the additional parameters 7 = (7j)j>i, as introduced in [29], have been 
found to be very useful from both a theoretical and a practical point of view. They allow 
us to associate more importance with some variables than with others, which ties in nicely 
with concepts such as effective dimension, see, e.g., [2]. This has been used to explain the 
success of quasi-Monte Carlo rules in finance. 

We now give the structure of the paper. In Section [2l we recall the definition of 
polynomial lattice rules, the scrambling algorithm and define the function space studied 
in this paper. The variance of estimators based on scrambled polynomial lattice rules is 
studied in Section O Next, in Sections S] and we construct polynomial lattice rules 
for which the variance of the associated estimator converges at a rate of iV-(i+2")+'5, for 
all S > 0. The constructions are based on a component-by-component approach and the 
Korobov construction respectively. In Section [6l we define a large class of randomized 
algorithms, which includes adaptive ones, and consequently establish that the variance of 
any estimator based on an algorithm from this class converges at most a rate of A^~(^+2") 
for the function space under consideration in this paper. Hence our constructions are 
almost optimal for the class of algorithms defined in Section [61 In Section [3, we study the 
implementation of the component-by-component approach, in particular, we show how 
to reduce the computational effort associated with it. This implementation is made use 
of in Section [HI where we compare the performance of scrambled polynomial lattice rules 
constructed in Section [H to the performance of scrambled digital nets. 

2 Preliminaries 

In this section, we define polynomial lattice rules, recall the scrambling algorithm and 
introduce the function space under consideration in this paper. 
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2.1 Polynomial Lattice Rules 

Polynomial lattice rules were introduced in [19], see also O [71 [18]. We fix a prime b and 
denote by the finite field containing b elements and by Z,h{{x~^)) the field of formal 
Laurent series over Z^. Elements of Zb((x~^)) are formal Laurent series, 

oo 
l=w 

where w is an arbitrary integer and all ti E l^b- The field Zb((a;~^)) contains the field of 
rational functions over Zf, as a subfield. Finally, the set of polynomials over Zf, is denoted 
by Zb[x]. For an integer m, we denote by Vm the map from Zb{{x'^)) to [0, 1) defined by 

' oo \ m 

Y^tix~'\= tib-'. 

^l=w / i=max(l,iii) 

The following definition of polynomial lattice rules stems from [in] , see also [3 [TH] . 

Definition 2.1 Let b be prime and m be an integer. For a given dimension s > 1, 
choose p{x) G Z,i,[x] with deg{p{x)) = m and . . . , qs{x) G Z(,[x]. For < h < b"' let 

h = hQ + hib + - ■ ■ + hm~ib™'^^ be the b-adic expansion of h. With each such h we associate 
the polynomial 

m—l 



h{x) = hrX^ G Zf,[x] 



r=0 

Then Sp.m{q), where q = (gi, . . . , g^), is the point set consisting of the b^ points 

for < h < b"^. A quasi-Monte Carlo rule using the point set Sp^raio) is called a polyno- 
mial lattice rule. 

We remark that polynomial lattice point sets are also digital nets, see [TJ [TTJ [18] . 

For the remainder of the paper, we use the following notation: We write h for vectors 
over Zfe and h for vectors over Z or M. Polynomials over Zb are denoted by h{x) and 
vectors of polynomials by h{x). Furthermore, given an integer h with 6-adic expansion 
h = Yl'^o hrb^, we denote the associated polynomial by 



h{x) = h 



T 

r=0 



T 

X 
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For arbitrary h{x) = {hi{x), . . . , hs{x)) G Zb[xY and q{x) = {qi{x), . . . , qs{x)) G 
we define the "inner product" 

s 

h{x) ■ q{x) = ^ h,jq.j{x) G 'Li\x] 
i=i 

and we write q{x) = (mod p{x)) if p{x) divides q{x) in Z,i,[x]. 

Finally, we introduce the dual lattice which plays an important role in numerical 
integration, see [5],[7j, which requires us to introduce the following function: For a non- 
negative integer k with 6-adic expansion k = ko + kib + . . . we write tTm{k) = ko + kib + 
■ ■ ■ + km-ib"^~^ and thus the associated polynomial 

trrnik){x) = ko + kix + . . . km-ix"''^ G Zfefo;] 

has degree < m. For a vector k G Ng, trm(fc) is defined componentwise. 

Definition 2.2 Let q{x) = (qi{x), . . . , qs{x)) G 'Z'b[xY , then the dual polynomial lattice of 
(q) is given by 

^ = %{q) 

tTmiki){x) H h timik2){x)q2 H h tTm{ks){x)qs{x) = (mod pix))}. 

2.2 The Scrambling Algorithm 

The scrambling algorithm is a randomization algorithm which was introduced by Owen, 
see [23] and also P [lOl [SI [251 [301 [3ll [33] . 

We now describe the scrambling algorithm using a generic point x G [0, 1)*, where 
X = (xi, . . . , Xg) and 

Then the scrambled point shall be denoted by j/ G [0, 1)"*, where y = {yi, . . . ,ys), 

The permutation applied to ^^y, j = 1, . . . , s depends on fc, ioi 1 < k < I. In particular, 
Vj,i = 7rj(^j,i), Vj,2 = -^j,^,,i{^j,2), Vj,3 = and in general 
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where tij and ^-i^ k > 2 are random permutations of {0, 1, . . . , 6 — 1}. We assume 

that permutations with different indices are mutually independent. Also, if we apply the 
scrambling algorithm to x to obtain y, then y is uniformly distributed in [0, 1)*, see [23| 
Proposition 2]. Finally, it was shown in |23] that the scrambling algorithm preserves the 
(t,m, s)-net property with probability 1, i.e applying the scrambling algorithm to the 
points of a (t, m, s)-net results in a {t,m, s)-net with probability 1. 

2.3 A Weighted Walsh Function Space based on Variance 

In this section, we introduce the function space under consideration in this paper. In 
particular, we consider weighted spaces and for this purpose, we introduce a sequence of 
positive, non-increasing weights 7 = (7^)°?,^ The purpose of the weights is to model the 
importance of the different variables and we point out that the idea stems from [29]. For 
s G N, let [s] = {1, . . . , s} and for u C [s] let 7^ := H^gu^i weight associated with 

the projection onto coordinates whose index is contained in u. 

Walsh functions have been an important tool in the analysis of digital nets; in |14j . 
Walsh functions were used for the first time to analyze nets and the connection between 
numerical integration using digital nets and Walsh functions was made in [6], see also [7]. 

We now briefly recall the definition of Walsh functions. Let Nq denote the set of 
nonnegative and N the set of positive integers. Each A; G N has a unique fe-adic repre- 
sentation k = Yli=o with digits fcj G {0, . . . , 6 — 1} for < z < a, where ka 7^ 0. For 
k = we have a = and ko = 0. Similarly, each x G [0, 1) has a 6-adic representation 
X = Yl'i^i with digits G {0, . . . , 6 — 1} for i > 1. This representation is unique in 
the sense that infinitely many of the must differ from h — 1. We define the fcth Walsh 
function in base 6, waU : [0, 1) — C by 

walfc(x) := exp(27ri(^i/co H h ^a+ika)/b) . 

For dimension s > 2 and vectors k = {ki, . . . ,ks) G Ng and x = (xi, . . . , Xg) G [0, 1)'^ we 
define walfc : [0, 1)" ^ C by 

s 

walk{x) := JJwaU^(Xj) . 

i=i 

Studying integration errors resulting from the approximation of an integral based on a 
digital net or a polynomial lattice rule, it is useful to consider the Walsh series of the 
integrand /. For / G i^2([0, 1]^), the Walsh series of / is given by 



fix) ~ /(fc)walfc(a;) , 

fceNf, 



(6) 
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where the Walsh coefficients f{k) are given by 



f{x)wa\kix)dx . 



[0,1] 



We do not necessarily have equality in Equation (|6]), however, the completeness of the 
Walsh function system {waU : k G Nq} (see for instance [71 Appendix A]), implies that 

Var[/]= J2 

ke%\{0} 

where Var[/] = L^.,{f{x) - f^dx, and where / = L^,, f{x)dx. 



Let c^(z,,o) 



'[0,1] 

^feei(u„o) 1/(^)1'' ^^^^^ 



L{Kfi) = {fc G Nq : 6'' ^ < fej < 6'" for z G u and /Cj = for z G [s] \ u} . 



(7) 



Further let = X]j=i h for Z = (/i, . . . , Is). For < a < 1 we define a weighted norm 
for functions / G L2([0, 1]"*) by 



max 7„ 



1/2 g,,p 6"l'"IV(,^.,o)(/). 



For < a < 1 define a space Va^s,'y ^ -^2([0, 1]**) consisting of all functions / for which 
II /II a < oo. (One could of course use some £p norm instead of the supremum-norm to 
define || ■ ||a and the function space, but these do not yield a quality criterion of the form 
([5]) which can be used for our construction, see (ITT]) and Lemma [3. II below.) 

The following observation stems from [TJ Section 13.5]: For a subinterval J = Y[i=i[^i^ Hi 
with < Xj < ?/j < 1 and a function / : [0, 1)'^ — > M, let the function A(/, J) denote the 
alternating sum of / at the vertices of J where adjacent vertices have opposite signs. 
(Hence for / = HjLi fi we have A(/, J) = ]Yi=i{fi{xi) - fi{yi)).) 

We define the generalized variation in the sense of Vitali of order < a < 1 by 



K^^)(/) = sup(^Vol(J) 



A(/, J) 



Vol(J)" 



1/2 



where the supremum is extended over all partitions V of [0, 1]'* into subintervals and 
Vol( J) denotes the volume of the subinterval J. 

For a = 1 and if the partial derivatives of / are continuous on [0, 1]** we also have the 
formula 

2 \ 1/2 



[0,1]- 



dxi ■ ■ ■ dxs 



dx 
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Until now we did not take projections to lower- dimensional faces into account. 
For 7^ it C [s], let Va^^\fu]u) be the generalized Vitali variation with coefficient 
< a < 1 of the Inl-dimensional function 



For M = we have = /^^ f{x)dx[s] and we define Va^^\ffD; 0) = I/0I. Then 



\uC[s] J 

is called the generalized Hardy and Krause variation of / on [0, 1]'*. 

A function / for which Vq(/) < 00 is said to be of finite variation of order a. 
The following result is from jl] and [3 Section 13.5]. 

Corollary 2.1 Let b > 2 be a natural number and let f G L2([0, 1]*) have bounded 
variation Va{f) < 00 0/ order < a < 1. Then 

ll/IU < max f ||/|U,70-\ VM) ,2^ax 7-^/'(& " lY^'^'^^^A . 

Hence every function / G ([0, 1]*) which has bounded variation of order < a < 1 
is in Va,s,-y For the extreme case a = one obtains Vo,s,7 = -^2([0, 1]'^), but this case is 
not included in our investigations since the criterion ([5]) is infinite in this case, see (fTTil 
and Lemma [3.11 below. 



3 The Variance of Estimators based on Scrambled 
Polynomial Lattice Rules 

In this section, we discuss the variance of the estimator 

A/) = ^E/W' (10) 

h=0 

where the points y^, . . . , ybm_i are obtained by applying the scrambling algorithm to a 
digital (t, m,s)-net over Zf,. 



10 



We use the following notation: For a non-negative integer k with 6-adic expansion 

k = ko + kib + . . . , 

we write k = (fco, fci, . . . ) , which is an infinite-dimensional vector, and we use 

tTm{k) = {ko, ki,..., fcm-l)^ • 

We now introduce the integration problem studied in this paper, in particular, we are 
interested in the worst-case variance of multivariate integration in V^^g^^ using a scrambled 
quasi- Monte Carlo rule Qb^,s'- 

Var(Q6m^s, Va,s,7) = sup Var[/(/, Qfem,^)] , 
/eK,,,,^,||/IU<i 

where /(/, Qb"^,s) denotes the estimator based on the point set obtained by applying 
the scrambling algorithm to Qb'^,s- We denote the quasi- Monte Carlo rule based on a 
polynomial lattice rule Sp^miq) by Qb"\s{q) and the associated worst-case variance by 
Var(Q6m ,,(qr), ^,,^,7)- For k = ko + Kib -\ h Ka-ib"'^ G No let 



ifk = 



^ if A: > 0. 



For k = {ki, ...,ks)eWQ let r„,^(fc) = n5=i ra,^,{kj). 

The next corollary gives a bound on the quantity ya.T{Qbm^s{q),Va,s,-y)- 

Corollary 3.1 Let < a < 1, q G l^bl^Y be a generating vector for a polynomial lattice 
rule with modulus p, and VeiT{Qb^,s{(l),Va,s,'y) be defined as above. Then 

VaT{Qb^,s{q),Va,s,f) < ^ r2a+l,7(fc), 

fcG5^p{q)\{0} 

where S>p{q) is the dual polynomial lattice. 

Proof. The corollary follows from the following facts: For any / G L2([0, 1]"), let /(/) 
be given by Equation (|TOl) and {xq, . . . ,Xbm_i} be a digital (t, m, s)-net over with 
generating matrices Ci, . . . , over Zf,, then we have (see [31] or also [H Section 13.5]) 

v-[A/)]= E E %^l^(..o)n^(C„...,c.)|, 
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where &{Ci, . . . , C,) = jfc e : Cjtimih) + ■■■ + Cjtr„(^,) = o} and L(j^„o) is given 
in Equation ([7]). Furthermore, if the Ci, . . . ,Cs are the generating matrices of the point 
set Sp^rniq), then for any G Nq \ {0} we have 

Cjtr„.{ki) + ■■■ + Cjtimih) = ^ tr^(fc) ■ g = (mod p) , 
which was first estabhshed in [T8| Lemma 4.40]. Using ([8]) and 

^|u| 

E (b-l)\n\ E fc(2.+l)|^.,|j A^u.0)n^p(g)|= Yl ^2a+l,-,(fe), 

the result follows. □ 



We denote the bound in Corollary 13.11 by 

B{q,a,j):= ^ r2„+i,^(A;). (11) 

keS>p{q)\{0} 

This bound is almost the same as the square worst case error for integration in a certain 
Walsh space considered in [5], see in particular [H Lemma 4.1]. 

As in [5], B{q,a,j) can easily be computed and therefore be used as a quality cri- 
terion for polynomial lattice rules. We write logj, for the logarithm in base b and we set 

^2aLlog6 0j ^ Q_ 



B(q, a, 7) = -1 + ^ 5^ n ( 1 + I 



Lemma 3.1 Let B[q,a,^) be given by Equation f lTT]) . T/ien 

ft™-! s 

' ' '12) 

\ 6 - 1 y ' 

where for x E [0, 1) we set 

5 _ 1 _ 52a[logjxJ(^^2o+l _ ]^-) 

'^-^^^ = Kfe^^^^T^ • 

A detailed proof of this result can be found in [1]. 

In the next remark, we show that if we construct a polynomial lattice rule which 
achieves optimal convergence rates for functions in V^^g^^ for some given < a < 1, then 
this polynomial lattice rule also achieves optimal convergence rates for functions in Va'^s,-/' 
where a < a' < 1. This means that the polynomial lattice rule constructed to achieve 
optimal convergence rates for functions of smoothness a adjusts itself to the optimal rate 
of convergence, as long as the smoothness a' of the function under consideration satisfies 
a' > a. 
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Remark 3.1 Assume that for a fixed a, < a < 1, we have constructed a polynomial 
lattice rule Sp^mio) such that 

B{q,a,j)<Cs,a,^,sN^^'^^''^^\ (13) 

for all 6 > 0, where Cs,a,-y,5 is permitted to depend on s,a,7 and 6. We point out that 
explicit constructions of polynomial lattice rules satisfying Equation f|T3|l are given in 
Sections\^and\B It follows immediately from Jensen's inequality, that 

l + 2a' . l + 2a' 

for a < a' < 1. Making use of Assumption f|T3l) . we conclude that 

for all 6 > 0. In particular, this observation motivates the construction of polynomial 
lattice rules for which a < 1, as the resulting point sets still achieve optimal convergence 
rates for functions of bounded variation of order a' , where a < a' < 1. 

4 Component-By-Component Construction of Poly- 
nomial Lattice Rules 

In this section, we show how to construct a polynomial lattice rule using a component- 
by-component approach so that the bound given in Equation ffTTl) converges at a rate 
of jv~^~^"+'', for any 6 > 0. We remark that in [71 Theorem 13.24], the corresponding 
result for digital nets was presented. A component-by-component (CBC) approach was 
first considered in [28] in the context of constructing lattice rules. Subsequently, the CBC 
algorithm has been applied to the construction of polynomial lattice rules in [5]. 

We use Rb,m to denote the set of all non-zero polynomials in with degree at most 
m — 1, i.e. 

Rb,m ■= {q £ '^b[x] '■ deg{q) < m and 7^ 0} . 

It is clear that |-Rfe,m| = &™ — 1 and furthermore it follows from the construction principle 
described in Subsection 12.11 that the polynomials qj can be restricted to Rb,m- Algorithm 
[1] gives the CBC algorithm. 

The next theorem shows that Algorithm [1] indeed constructs a G -R^^ so that 
B{{ql, . . . , g^), a, 7) converges at a rate of ]\[-'^-'^"'+^ ^ for any 6 > 0. 
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Algorithm 1 CBC algorithm 

Require: b a prime, s,m £ N and weights 7 = (7j)j>i- 
1: Choose an irreducible polynomial p G ^^[x], with deg(p) = m. 
2: Set gi = 1. 
3: for d = 2 to s do 

4: find G -R;,,m by minimizing -B((gi, . . . , g^), a, 7) as a function of g^. 
5: end for 

6: return q = (gi, . . .,qs). 



Theorem 4.1 Let b be prime and p G Zfe[a;] 6e irreducible, with deg{p) = m > 1. Suppose 
{ql, . . . , g*) G -R^„^ zs constructed using AlgorithmUl Then for all d = 1, . . . ,s we have 

B{{ql . . . , g^), a, 7) < .^,„_\.VA n [1 + ^M.a] , 



V ' j=i 



/or a// < A < 1 where 



( 1 

a,a,A = max (^ ^^2a_l)A ' ^2aA _ J " (^4) 

A proof of this result can be obtained by making a few modifications to the proof 
of [H Theorem 4.4], which is included in [1]. The additional term (6^" — l)^'^ in the 
definition of Cb,a,x arises from the one-dimensional case for which we have (we assume 

k = Ko + Kib-\ h Ka-lb°-~^) 

, oo 

k=l,b"^\k 
, oo 

l=m+l 
1 7l 



^(2a+l)m _ I 



for all < A < 1. The induction with respect to the dimension can be carried out as 
in the proof of O Theorem 4.4]. 

The next result discusses the tractability of Algorithm [T] 
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Corollary 4.1 Let b be prime, p G 1^b[x] be irreducible with deg{p) = m > 1 and N = b"" 
Suppose E Rl^ is constructed using AlgorithmUi Then we have the following: 



1. 

where 
2. Assume 



B{ql, a, 7) < €^^^^,5 - l)-(2°+i)+^ , for all < 5 < 2a, 

2q+1-(5 



1 

J- + Ij <-^f),a,(2a+l-5)-i 



Then Cs^a,-y,5 < Coo,a,7,5 < oo and we have 

B{q*s, a, 7) < c^^^^^^N - i)-(2°+i)+5 ^ y^r a// < 5 < 2a. 
Thus the bound B{q*,a,^) is bounded independently of the dimension. 
3. Under the assumption 

A := lim sup — < cxo, 



A+ri 

we obtain Cg^a,-y,2a < c^{b — l)^"s''^"-i and therefore 

A+ri 

Biq:,an)<5rjS^^iN-iy' 



for all rj > 0, where the constant only depends on rj. Thus the bound B{q*,a,^) 
satisfies a bound which depends only polynomially on the dimension. 

The proof is similar to the proof of [5], Corollary 4.5] and can be found in [1]. 



5 Construction of Korobov Polynomial Lattice Rules 

In this section, we construct Korobov polynomial lattice rules. The ideas underlying 
this algorithm stem from the construction of lattice rules, see [12]. We remark that 
the construction of Korobov polynomial lattice rules has been examined in [5], see also 
[13]. We denote the generating vector for the Korobov polynomial lattice rule by ip^q) = 
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Algorithm 2 Korobov algorithm 

Require: b a prime, s,m G N and weights 7 = (7j)j>i- 
1: Choose an irreducible polynomial p G ^^[x], with deg(p) = m. 
2: Find g* G -Rfe,m by minimizing ^(^/^(g), a, 7). 



(1, g, . . . , g*~"^) (mod p). As in SectionlU we work with the bound B{ilj{q), a, 7) and now 
state the algorithm showing how to construct Korobov polynomial lattice rules. 

We obtain the following bound for B{ip{q*),a,j), where q* is constructed using Algo- 
rithm[2l The proof of the following theorem can be obtained by making a few modifications 
to the proof of [5l Theorem 4.7], which are presented in [1]. 

Theorem 5.1 Let b be prime, s > 2 and letp G be irreducible with deg{p) = m > 1. 
A minimizer q* obtained from Algorithmic satisfies 



for all ^ < X <1, where Cb,a,x > is given by (fT4l) . 

We point out that the bounds in Theorems 14.11 and 15.11 only differ by the additional 
factor s^^^. We remark that the same observation was made in [S] and is also known 
from the lattice rule case. This leads to the conclusion that the Korobov construction is 
inferior to the component-by-component construction. 

In the next corollary, we discuss the tractability of Algorithm [2l 

Corollary 5.1 Let b be prime, s > 2, p E Zf,[x] be irreducible with deg{p) = m > 1 
and N = . Suppose q* G Rb,m is constructed using Algorithmic Then we have the 
following: 



B{ij{q*),a,-f) < 



n (1 + 7,'a,.,A) 



s 



1/A 



- l)i/A 



1. 



B{^{q*),a,j)<c. 



^2a+l-5(^ _ ^yi2a+l)+5 ^ all < 6 < 2a, 



where 



c 




2. Under the assumption 



A : = lim sup, 



< oo 



log s 
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we obtain 

A+r, 

and therefore 

for all rj > 0, where the constant only depends on rj. Thus the hound B{tlj{q*), a, 7) 
satisfies a bound which depends only polynomially on the dimension. 

The proof is again similar to the proof of Corollary [SI Corollary 4.8] and can be found 
in p. 



6 A Lower Bound on the Worst-Case Variance 

In this section, we produce a lower bound on the worst-case variance discussed in Sec- 
tion [3l As we rely on pOl Section 2.2.4, Proposition 1] to establish the result, the class of 
algorithms to which our result applies is the same as the class considered there. We now 
recall the definition of this class. Following [Sni Section 1.1], we use the notation 

S{f) = [ f{x)dx, 

for / e V^,s,7 and consider approximating S : Va^s,-y ^ using a mapping S : V^,s,-y — > M. 
As in [20I Section 1.1], we assume that in general, the function / G V^,^,^ is not known, 
but we have some information on / available, which is denoted by L, where L : K«,s,-y — )■ H 
and an approximation S : Va,s,-y K only uses the information L if it can be written as 
follows S = (foL, where M is a an arbitrary mapping, referred to as an (idealized) 

algorithm in [20] . In particular, we allow our approximation nodes to be chosen adaptively 
and define the following information operator: 

In = {L: ^ = (/(ai), fia^UM]), /(ajv[/(ai), . . . , /(a^-i)])) , 

where ai G [0, 1]" and : R*~^ [0, 1]' for i = 2, . . . , s} 

and we can now introduce the class of all approximations considered in this section: 

An = : Va,s,'y ^R\S = <^oL with : ^ M and L e /jv} . 

We remark that non-adaptive algorithms are of course included in An, consider S = ipoL, 
where L{f) = (/(ai), . . . /(aAr)). Now, following [Sni Section 2.1], we can define the 
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randomized algorithms considered in this paper, referred to as generahzed Monte Carlo 
methods in [20]: A random variable Q = {Q{uj))i^^n is called a randomized algorithm 
in An if {fl,B,fi) is a probability space and Q^u) G Ajy for all u E fl. The set of 
all randomized algorithms is denoted by *C{Ai\f), hence randomly scrambled nets (and 
therefore polynomial lattice rules) are also included in this set. We now present the lower 
bound on the worst-case variance, which applies to all randomized algorithms in *C{Aj^). 

Theorem 6.1 Let *C{An), Va,s,-y be defined as above. Then 



Proof. We remark that this proof follows along the lines of the proof of [H Theorem 10]. 
We only consider s = 1, since integration in V^^i^^-^ is no harder than integration in Va^s,-f 
with s > 1, as the one-dimensional space Va^i,^-^ can be identified with the subspace of 
Va^s,-^ consisting of functions depending only on the first variable. We let N be any given 
natural number and choose an integer m such that 



inf sup Var(J(/,g)) > CN' 



for some constant C independent of N where 




b'- 



,m— 1 



<2N <b'' 



m 



We define basic intervals 




and let ga{x) = 1b^.„(x) be the characteristic function of Bm,a- Then 




We now define 



ft"*-! 




a=0 
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where C,a G {1, —1} and bound (jf{g). Using Plancharel's identity we obtain that for any 
/ > we have 



crfig) < '^(rf,{g) = / g'^{x)dx = ^ ^^^^ / ga{x)g^{x)dx = -r7^^C = ^ 

l'=0 a,c=0 "^0 ^ a=0 

Further, for k > b"^ we have 



/■i " /"i " /"(.a+ij/b"' 

/ 5((x)walfc(x)rfx = ga{x)walk{x)dx = ^ walk{x)dx = 

-^0 a=0 a=0 "^«/^'" 



since j'^^l^^^ wslk{x)dx = for A; > 6"* and hence for / > m we have 

^K9)= E 1^(^)1' = 0- 

We set fa = 7i&~""*5'a for a = 0, 1, . . . , 6™" — 1. These fa have disjoint support and 

fa{x)dx > 7ir("+l)"^. 

[0,1] 

Set 

ft"*-! 



a=0 



then we get crfif) < 7i& for < / < m and cr;^(/) = for / > m. Hence 
||/|U = 7r^up6"'a,(/) <7fi sup 6"S&-""^<1 

iGN l<Z<m 

and the resuh follows now from [201 Section 2.2.4, Proposition l(ii)]. □ 

Remark 6.1 For a large class of randomized algorithms, including adaptive ones, we 
have shown that the worst-case variance in the Walsh function space Va^s,-/ behaves like 
j\j-{i+2a) ^ /n Sections \^ and we presented two algorithms which achieve worst-case 
variances of order N~^^~^'^°'^~^^ , for all S > 0, and are hence almost optimal for the class 
of algorithms *C{Aisf). 
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7 Implementation of the Component-By-Component 
Algorithm 

In this section, we show how to implement the CBC algorithm from Section HI Our 
approach is based on [22], but we simplify the algorithm using ideas from [3]. Using ideas 
from [211 [22], we obtain, for d > 2, 



h=0 j=i 
d 



b \ 1 ^""^ / b 

^ h=l ^ 



6-1 



Xh.d, 



where 



i=i ^ ^ 



Let u (^j = (t)o.{xh,d), where h and denote the polynomials associated with h and 

Qd and p denotes the polynomial p = p{x) G ^^[x]. Following [21], we now introduce the 
following matrix 

hq^ 



00 



p 



q=l,...,b"^-l 
/i=l,...,b'"-l 



(16) 



i.e. rows are indexed by q and columns by h. 

Let Pd-i = (Pd-i(l)) • • • ^Pd-iiP^ ~ ■ Following [21], we have an update rule for 
Pd given by 

Pd = diag (^(^l(6-_i)x(6--i) + -g-—[^d^^ 

where diag(a;) denotes the diagonal matrix with the elements of x on its diagonal and 
zero elsewhere and where we use Vj to denote a selection vector with 1 in position j and 
elsewhere. 

We now use the notation B^^i = (i?((q^_]^, 1), a, 7), . . . , 5((q^_;^, 6"^ — 1), a, 7))^. 
Then 



Bd-i 



L(b'"-l)xl 
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1 / b \ 



-{6'"-l)xl 



/l=l 



In the next lemma, we summarize an observation from [3]. Let 

Il{g) = [Ilk,i\k=i,...,b"'~i 
;=i,...,fe™-i 



where 



and 



where 



1 if k{x) = (mod p) 

otherwise 



l=l,...,b"^-l 



1 if = g \x) (mod 
otherwise, 



(17) 



be two permutation matrices, where g is a. primitive element which generates all elements 
of (Zf,[x]/p)* = {g^,g^, ■ ■ ■ jf?^""^}; such an element g is known to exist since the multi- 
plicative group of every finite field is cyclic. Let = deg{g'' (mod p)), k = 0,1, . . . , b"^ — 2, 
and set 

A3 = [6^°*^-^- (-°'i''"-i)]i=i,...,b--i (19) 

j=i'...'fe'"-i 

and note that is a circulant matrix, which allows us to use Fast Fourier Transforms 
(FFTs) as in [211 [22] . We now state the lemma. 

Lemma 7.1 Letp be an irreducible polynomial, let g be a primitive element of {Zh[x]/p)* , 
and let Il{g), Il{g~^), As and flp be defined as above. Then 



b-1 



Qp - l(fe'»-l)x(fe'»-l)^^^2a _ n - h(h2a _ ^ 



-2am 



'6(62^-1) b{b^ 



Il{g)AsU{g 



-1\T 
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Proof. It follows from the definition of (paix) 

( f^^W - (52a+l _ ^)^-2aao.h., 

V " V7;; " b{b^- - 1) 6(62- _ 1) ' 

where ao,/i,g denotes the smallest integer a so that ^h,q,a 7^ 0, and where 

hq\ ^h,q,l , ^h,q,2 , 



g / 6 6^ 
Hence 

5 _ 1 52a+l _ 1 

where 

/i=l,...,b'"-l 



Now assume that for w G Zb[x] we have 

+ U2,wX~'^ + (20) 



p{x) 

where Uj^y^ ^ ^6- Then 

) = Wi,7^^6-' + u, ji^b-^ + ■■■ + u^jJo-^ , 

hence, ao^h,q is the smallest integer a so that w^j^g 7^ 0, /i, g G (note that p J(h,q). 

The matrix given by 

A2 = Il'ig)A,Uig~') 
is circulant. Indeed it can be checked that 

^2 = r6-2"%,.-,«'],=i,„„,™_i, (21) 

J=l,...,b'"-1 

where and in Equation (12T|) denote the integers associated with the polynomials 
(mod p) and (mod p). We let a^^g-j^gi = Vi^j and note that = r^/ for k = k' 
(mod 6™ — 1), as g^"""^ = 1 (mod p), hence 

j=l,...,6'"-l 
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The matrix A2 is circulant and is the smallest integer r such that Uj.^gk 7^ 0, which 
implies using Equation ( l20l) that 

deg(5f^ (mod p)) = deg(p) - , 

= m — deg{g^ (mod p)) . 
tk = deg{g'' (mod p)) , 



and consequently 
Now denoting 
we get 



12 — w 7-13 , 

where A3 is given by Equation f|T9|) and the result follows. □ 

Note that if the polynomial p in the lemma above is primitive, then one can choose 
the primitive element g{x) = x. In Algorithm |3] we show how to implement the CBC 
algorithm from Section HI Several remarks regarding Algorithm |3] are in order. 

Algorithm 3 Fast CBC algorithm 
Require: b a prime, s, m G N and weights 7 = (7j)j>i. 
1: Choose a primitive polynomial p G ^^[x], with deg(p) = m, and choose g{x) = x. 

Mo l(6'"-lxl)- 

for d = 1 to s do 

Wd = argmin^gRj,,^^ Bdiw) . 



= diag (^llx(fe'"-l) \l + ^^2a _ 1) J - ^'^ (5_l)(52a_l) -^3(^^d, J Z^rf^l 



7: end for 

8: return q = (gi, . . .,qs). 



Remark 7.1 ^45 in f21\ l2^ . we search for the minimum in the permuted space, hence we 
minimize Bd = 11^ {g)Bd^i. However, as in f2l\ \2^ . the component Zd can he found by 
mapping back Wd using Tl{g). 
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Remark 7.2 We have fx^ = 11^ {g~^)p^ and consequently update fi^ using ^^-i- Hence 
we do not need to permute hack and forth, hut can complete the algorithm in the permuted 
space. 

The next corollary gives information on the computational complexity of Algorithm 
[31 We use 



a 



(22) 



to denote the vector generating the circulant matrix A3 in Lemma 17.11 and Algorithm |3l 

Corollary 7.1 Assume that the vector a in Equation fl22l) has heen precomputed and 
stored using 0{b"^) memory. Then Algorithmic can he completed in time 0{sb'^m) and 
memory 0{b"^). 



For a proof, see [2T1 [22] or also [71 Section 10.3]. 



8 Numerical Experiments 

In this section, we numerically investigate the performance of the CBC algorithm pre- 
sented in Section [H we rely on Section [71 for the implementation of the algorithm. In 
Tables m - [3l we present values of B{q,a,^) for different choices of a and 7, where q is 
constructed using Algorithm [31 

We compare the performance of the CBC algorithm to the performance of digital 
nets. As was done with scrambled polynomial lattice rules in Section [31 we can study the 
variance of the estimator /(/) given in Equation (ITOl) . consider the worst-case variance of 
multivariate integration in Va,s,-y and bound this variance as follows: 

Var((5f,m,s(Ci, ...,Cs), Va,s,f) < ^ (k), (23) 

fce5^(Ci,...,c,)\{o} 

where Ci, . . . ,Cs are the generating matrices of the digital net under consideration and 
^{Ci, . . . ,Cs) is its dual space. We denote the bound f[23l) by i?((Ci, . . . , C^), a, 7), 
and remark that B{(Ci, . . . , Cs), a, 7) can also be computed using Equation f[T2l) . where 
{xhjj^^Q is the digital net generated by Ci, . . . , Cg- 

Consequently, we compare the values ofB{q, a, 7) to the values of B{{Ci, . . . , Cg), a, 7) 
in Tables [H - [31 in each cell, the top number corresponds to the CBC construction and 
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the bottom one to the digital net. We choose the following digital nets: For s = 1, we 
simply choose equidistributed points, Xh = h = 0,...,b"^ — 1, for s = 5, we use 
Pirsic's implementation of Niederreiter-Xing points, [26], and for s = 50 and s = 100, we 
use Sobol points as constructed in [11]; we point out that for the CBC construction, we 
choose b = 2 and likewise, the digital nets under consideration are digital nets over Z2. 

We derive the following conclusions from the tables: For s = 1, as expected, we 
obtain the optimal rate of convergence, 2"'^^°"'"^)'", and observe the same values for the 
CBC construction as for the digital nets. Regarding the case s = 5, the values are 
comparable, however, the Niederreiter-Xing construction seems to be slightly better than 
the CBC construction for the examples considered. Finally, for s = 50 and s = 100, 
the performances of the two methods are again comparable, however, this time, the CBC 
construction seems to outperform the digital nets. 
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a = 


= 0.5 






a 


= 1 






m = 


s = 1 


s = 5 


s = 50 


s = 100 


s = 1 


s = 5 


s = 


50 


s = 100 


A 


3.91e-03 


1.46e+00 


7.04e+13 


7.92e+28 


8.14e-05 


4.37e-02 


l.lOe- 


f05 


1.95e+ll 


3.91e-03 


1.48e+00 


7.04e+13 


7.92e+28 


8.14e-05 


4.90e-02 


l.lOe+05 


1.95e+ll 




9.77e-04 


6.16e-01 


3.52e+13 


3.96e+28 


1.02e-05 


1.09e-02 


5.52e- 


f04 


9.74e+10 


o 


9.77e-04 


6.34e-01 


3.52e+13 


3.96e+28 


1.02e-05 


1.32e-02 


5.52e- 


f04 


9.74e+10 


u 


2.44e-04 


2.66e-01 


1.76e+13 


1.98e+28 


1.27e-06 


3.45e-03 


2.76e- 


f04 


4.87e+10 


2.44e-04 


2.61e-01 


1.76e+13 


1.98e+28 


1.27e-06 


3.17e-03 


2.76e- 


f04 


4.87e+10 


7 


6.10e-05 


1.08e-01 


8.80e+12 


9.90e+27 


1.59e-07 


9.05e-04 


1.38e- 


f04 


2.44e+10 


6.10e-05 


1.04e-01 


8.80e+12 


9.90e+27 


1.59e-07 


7.19e-04 


1.38e- 


f04 


2.44e+10 


c 


1.53e-05 


4.24e-02 


4.40e+12 


4.95e+27 


1.99e-08 


2.36e-04 


6.90e- 


f03 


1.22e+10 


o 


1.53e-05 


3.93e-02 


4.40e+12 


4.95e+27 


1.99e-08 


1.48e-04 


6.90e- 


f03 


1.22e+10 


Q 


3.81e-06 


1.74e-02 


2.20e+12 


2.48e+27 


2.48e-09 


6.10e-05 


3.45e- 


f03 


6.09e+09 


3.81e-06 


1.44e-02 


2.20e+12 


2.48e+27 


2.48e-09 


2.86e-05 


3.45e+03 


6.09e+09 


1 n 

iU 


9.54e-07 


6.41e-03 


l.lOe+12 


1.24e+27 


3.10e-10 


1.29e-05 


1.72e- 


f03 


3.04e+09 


9.54e-07 


5.21e-03 


l.lOe+12 


1.24e+27 


3.10e-10 


5.56e-06 


1.72e- 


f03 


3.04e+09 


1 1 


2.38e-07 


2.29e-03 


5.50e+ll 


6.19e+26 


3.88e-ll 


2.56e-06 


8.62e- 


f02 


1.52e+09 


2.38e-07 


1.82e-03 


5.50e+ll 


6.19e+26 


3.88e-ll 


l.Ole-06 


8.62e- 


f02 


1.52e+09 


12 


5.96e-08 


8.39e-04 


2.75e+ll 


3.09e+26 


4.85e-12 


5.03e-07 


4.31e- 


f02 


7.61e+08 


5.96e-08 


6.17e-04 


2.75e+ll 


3.09e+26 


4.85e-12 


1.78e-07 


4.31e+02 


7.61e+08 


13 


1.49e-08 


3.09e-04 


1.37e+ll 


1.55e+26 


6.06e-13 


1.05e-07 


2.15e- 


f02 


3.81e+08 


1.49e-08 


2.06e-04 


1.37e+ll 


1.55e+26 


6.06e-13 


3.07e-08 


2.16e- 


f02 


3.81e+08 


14 


3.73e-09 


1.12e-04 


6.87e+10 


7.74e+25 


7.58e-14 


2.56e-08 


1.08e- 


f02 


1.90e+08 


3.73e-09 


6.76e-05 


6.87e+10 


7.74e+25 


7.57e-14 


5.17e-09 


1.08e- 


f02 


1.90e+08 


15 


9.31e-10 


3.66e-05 


3.44e+10 


3.87e+25 


9.27e-15 


4.98e-09 


5.38e- 


fOl 


9.52e+07 


9.31e-10 


2.18e-05 


3.44e+10 


3.87e+25 


9.33e-15 


8.54e-10 


5.39e- 


fOl 


9.52e+07 


16 


2.33e-10 


1.29e-05 


1.72e+10 


1.93e+25 


1.22e-15 


8.92e-10 


2.69e+01 


4.76e+07 


2.33e-10 


6.94e-06 


1.72e+10 


1.93e+25 


l.lle-15 


1.38e-10 


2.69e+01 


4.76e+07 



Table 1: Values of B{q,a,^) and -B((Ci, . . . , C^), a, 7) for 7^ = 1, j = l,...,s and q 
constructed using the CBC algorithm; the top number gives the value of B{q,a,^), the 
bottom the value of B{{Ci, . . . , Cs), a, 7). 
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a = 0.5 


a = 1 


m = 


s = 1 


s = 5 


s = 50 


s = 100 


s = 1 


s = 5 


s = 50 


s = 100 


A 


3.42e-03 
3.42e-03 


4.64e-01 
4.84e-01 


2.03e+01 
2.04e+01 


2.04e+01 
2.06e+01 


7.12e-05 
7.12e-05 


1.47e-02 
1.83e-02 


2.82e-01 
3.46e-01 


2.83e-01 
3.48e-01 


o 


8.54e-04 
8.54e-04 


1.87e-01 
1.95e-01 


9.99e+00 
l.Ole+01 


l.Ole+01 
l.Ole+01 


8.90e-06 
8.90e-06 


3.54e-03 
4.45e-03 


1.16e-01 
1.38e-01 


1.17e-01 
1.39e-01 


u 


2.14e-04 
2.14e-04 


7.75e-02 
7.46e-02 


4.91e+00 
4.94e+00 


4.95e+00 
4.98e+00 


l.lle-06 
l.lle-06 


1.04e-03 
9.29e-04 


4.78e-02 
5.30e-02 


4.82e-02 
5.34e-02 


7 


5.34e-05 
5.34e-05 


2.96e-02 
2.80e-02 


2.40e+00 
2.44e+00 


2.42e+00 
2.47e+00 


1.39e-07 
1.39e-07 


2.54e-04 
1.97e-04 


1.85e-02 
2.39e-02 


1.87e-02 
2.41e-02 


fi 

o 


1.34e-05 
1.34e-05 


1.17e-02 
l.Ole-02 


1.17e+00 
1.20e+00 


1.18e+00 
1.21e+00 


1.74e-08 
1.74e-08 


5.77e-05 
3.79e-05 


7.39e-03 
1.08e-02 


7.45e-03 
1.09e-02 


Q 


3.34e-06 
3.34e-06 


4.43e-03 
3.54e-03 


5.66e-01 
5.89e-01 


5.71e-01 
5.95e-01 


2.17e-09 
2.17e-09 


1.29e-05 
6.98e-06 


2.84e-03 
4.94e-03 


2.87e-03 
4.97e-03 


1 n 


8.34e-07 
8.34e-07 


1.56e-03 
1.22e-03 


2.72e-01 
2.88e-01 


2.75e-01 
2.90e-01 


2.72e-10 
2.72e-10 


3.03e-06 
1.28e-06 


1.08e-03 
2.24e-03 


1.09e-03 
2.26e-03 


1 1 


2.09e-07 
2.09e-07 


5.45e-04 
4.10e-04 


1.30e-01 
1.42e-01 


1.31e-01 
1.43e-01 


3.40e-ll 
3.40e-ll 


6.24e-07 
2.22e-07 


4.01e-04 
9.58e-04 


4.06e-04 
9.66e-04 


12 


5.22e-08 
5.22e-08 


1.93e-04 
1.35e-04 


6.20e-02 
6.73e-02 


6.26e-02 
6.80e-02 


4.24e-12 
4.24e-12 


1.16e-07 
3.79e-08 


1.49e-04 
3.59e-04 


1.51e-04 
3.64e-04 


13 


1.30e-08 
1.30e-08 


7.07e-05 
4.38e-05 


2.94e-02 
3.33e-02 


2.97e-02 
3.36e-02 


5.31e-13 
5.30e-13 


2.48e-08 
6.34e-09 


5.43e-05 
2.11e-04 


5.51e-05 
2.13e-04 


14 


3.26e-09 
3.26e-09 


2.27e-05 
1.40e-05 


1.39e-02 
1.58e-02 


1.40e-02 
1.60e-02 


6.62e-14 
6.62e-14 


4.56e-09 
1.04e-09 


1.99e-05 
8.23e-05 


2.02e-05 
8.31e-05 


15 


8.15e-10 
8.15e-10 


8.01e-06 
4.41e-06 


6.49e-03 
7.70e-03 


6.57e-03 
7.78e-03 


8.49e-15 
8.22e-15 


l.lOe-09 
1.67e-10 


7.15e-06 
4.44e-05 


7.26e-06 
4.48e-05 


16 


2.04e-10 
2.04e-10 


2.68e-06 
1.37e-06 


3.02e-03 
3.76e-03 


3.06e-03 
3.80e-03 


9.99e-16 
8.88e-16 


1.81e-10 
2.64e-ll 


2.57e-06 
2.14e-05 


2.61e-06 
2.15e-05 



Table 2: Values of B{q, a, 7) and B{{Ci, . . . , Cs), a, 7) for 7^ = 0.875-', j = 1, . . . , s and q 
constructed using the CBC algorithm; the top number gives the value of B{q,a,^), the 
bottom the value of B{{Ci, . . . , Cs), a, 7). 
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a = 


0.5 






a = 


= 1 




m = 


s = 1 


s = 5 


s = 50 


s = 100 


s = 1 


s = 5 


s = 50 


s = 100 


A 


3.91e-03 


2.75e-02 


4.75e-02 


4.90e-02 


8.14e-05 


7.68e-04 


1.80e-03 


1.88e-03 




3.91e-03 


3.20e-02 


5.97e-02 


6.17e-02 


8.14e-05 


1.27e-03 


4.40e-03 


4.62e-03 




9.77e-04 


8.98e-03 


1.78e-02 


1.84e-02 


1.02e-05 


1.48e-04 


4.88e-04 


5.20e-04 


o 


9.77e-04 


1.25e-02 


2.22e-02 


2.30e-02 


1.02e-05 


3.13e-04 


1.23e-03 


1.30e-03 


u 


2.44e-04 


2.95e-03 


6.29e-03 


6.56e-03 


1.27e-06 


3.37e-05 


1.20e-04 


1.31e-04 


2.44e-04 


3.20e-03 


7.23e-03 


7.66e-03 


1.27e-06 


3.62e-05 


2.47e-04 


2.89e-04 


7 


6.10e-05 


8.96e-04 


2.22e-03 


2.35e-03 


1.59e-07 


5.05e-06 


2.91e-05 


3.31e-05 


6.10e-05 


l.lle-03 


2.65e-03 


2.88e-03 


1.59e-07 


8.38e-06 


6.91e-05 


9.12e-05 


fi 


1.53e-05 


2.96e-04 


7.86e-04 


8.36e-04 


1.99e-08 


1.04e-06 


6.94e-06 


8.16e-06 


o 


1.53e-05 


3.44e-04 


l.OOe-03 


1.12e-03 


1.99e-08 


1.37e-06 


2.16e-05 


3.52e-05 


Q 


3.81e-06 


9.34e-05 


2.81e-04 


3.02e-04 


2.48e-09 


1.90e-07 


1.67e-06 


2.02e-06 


3.81e-06 


9.15e-05 


3.27e-04 


3.64e-04 


2.48e-09 


1.60e-07 


4.70e-06 


6.18e-06 


1 n 


9.54e-07 


2.78e-05 


9.54e-05 


1.03e-04 


3.10e-10 


4.07e-08 


4.20e-07 


5.14e-07 


9.54e-07 


2.69e-05 


1.19e-04 


1.32e-04 


3.10e-10 


2.49e-08 


1.83e-06 


2.37e-06 


1 1 


2.38e-07 


8.95e-06 


3.34e-05 


3.64e-05 


3.88e-ll 


6.34e-09 


9.59e-08 


1.21e-07 


2.38e-07 


8.25e-06 


4.10e-05 


4.71e-05 


3.88e-ll 


3.76e-09 


3.25e-07 


5.44e-07 


12 


5.96e-08 


2.68e-06 


1.18e-05 


1.30e-05 


4.85e-12 


1.30e-09 


2.36e-08 


3.03e-08 


5.96e-08 


2.54e-06 


1.46e-05 


1.68e-05 


4.85e-12 


6.34e-10 


1.16e-07 


1.64e-07 


13 


1.49e-08 


8.29e-07 


4.05e-06 


4.50e-06 


6.06e-13 


2.04e-10 


5.79e-09 


7.61e-09 


1.49e-08 


7.08e-07 


4.69e-06 


5.76e-06 


6.06e-13 


9.21e-ll 


2.45e-08 


5.16e-08 


14 


3.73e-09 


2.50e-07 


1.40e-06 


1.57e-06 


7.58e-14 


4.11e-ll 


1.40e-09 


1.88e-09 


3.73e-09 


1.97e-07 


1.60e-06 


1.98e-06 


7.57e-14 


1.29e-ll 


7.04e-09 


2.04e-08 


15 


9.31e-10 


7.70e-08 


4.91e-07 


5.54e-07 


9.27e-15 


6.15e-12 


3.45e-10 


4.79e-10 


9.31e-10 


5.59e-08 


5.50e-07 


7.09e-07 


9.33e-15 


1.74e-12 


2.52e-09 


4.52e-09 


16 


2.33e-10 


2.34e-08 


1.71e-07 


1.95e-07 


1.22e-15 


l.OOe-12 


8.51e-ll 


1.22e-10 


2.33e-10 


1.69e-08 


1.89e-07 


2.52e-07 


l.lle-15 


2.70e-13 


9.54e-10 


1.45e-09 



Table 3: Values of B{q, a, 7) and B{{Ci, . . . , Cs), a, 7) for 7^ = j = 1, . . . , s and q 
constructed using the CBC algorithm; the top number gives the value of B{q,a,^), the 
bottom the value of B{{Ci, . . . , Cg), a, 7). 
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